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SUMMARY 




Equations comprising a solution for the time spent by a 
satellite in umbra and penumbra are derived. Umbra and penumbra 
regions are assumed to be right circular cones, and shadowing bodies 
are assumed spherical. Pertinent constants are evaluated for shadowing 
by both the earth and the moon. The solution is applied in predicting 
shadow time distributions for the AIMP mission (establishment of a 
lunar satellite orbit) . 
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A INTRODUCTION 


la connection with the AIMP mission, a requirement was manifest 
to predict shadow time distributions for lunar satellite orbits. At the 
t ime of the study no lunar orbit shadow program was available to the Delta 
Project, and tabulated data available were based on solutions which assumed 
cylindrical shadow regions. 

While cylindrical shadow theories are sufficiently accurate for most studies 
involving low or medium altitude satellite, accurate determination for 
high altitudes should be based on the conical model of umbra and penumbra 
regions . 


The order of magnitude of the difference in the two solutions is 
indicated by Figure 1. This figure shows the maximum time in umbra and the 


of various radii. Times for both earth 
satellites and lunar satellites are shown.. 
Times determined by a cylindrical shadow 
theory would lie between the "umbra only" 
and "umbra and penumbra" curves. The. 
separation of these curves can be seen to be 
relatively large for high altitude orbits. 



Because of error magnitudes possible in the injection into earth-to- 
moon transfer orbits, lunar orbits achieved in the AIMP mission may have 
an extremely wide range of eccentricities and altitudes. For this reason, 
a conical shadow theory was regarded as desirable. 


A solution in which the umbra region is assumed conical, and the 
shadowing body is assumed spherical, is presented in Section B. Appendix C 
shows the applicability of this solution to the determination of time 
spent in the penumbra. The. iterative solution is relatively simple and 
requires as inputs the. size and shape of the satellite orbit and the angular 
position of the sun relative to the orbit plane. 


Several authors 
(References 1 - 5) . 
However, Fixler (Ref. 
which is reducible to 


have considered the shadow time determination problem 
Most solutions assume the cylindrical shadow model. 
4) determines, by a different derivation, a solution 
that of Section B. 
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B. ANALYSIS 


The geometry o.f the shadow time computation problem is shown in 
Figure 2. Here (A) is an edge view of the orbit plane viewed along a 
line perpendicular to the moon-sun line, and (B) is a projection showing the 
top view of the orbit plane. Scale is distorted, of course, for clarity. 

The lunar orbit plane intersects the shadow cone in a conic, generally an 
ellipse, hereafter called the "shadow ellipse" 


If p is the umbra half -cone angle, 

sin p = - R f = 0.004646303 . 

Mo 6 

The value of MS here is taken as the average value, 149.53 x 10 km 
^see Appendix A) . 

p = 0 . 2649 & 
y,= 90° - S - p 
%= 90° - i + P 

Then the semimajor axis of the shadow ellipse, a s , can be obtained from 
II* =. COS y, = sin ( $+ p ) 
a nd p- 

^ = cos y a = sm(g-p) 

so that 

„ = AlL = R«r ! + — ! " . 

5 2 2 L sin ( S + p) sin ( S - p) 

Note that since 


sih ( S + p) sin (S - pj = sin 2 S — sin a p , 


( 2 ) 


= R Sir'S cos p (3 

5 h • 2 C *2 

sin i — sin p 

The shadow cone touches the moon along a minor circle of "latitude" p 
above the great circle normal to the moon -sun line. This minor circle 
touches the shadow ellipse in the orbit plane at two points, P and Q. The 
inclination of the lunar orbit plane to the normal, great circle is 90° - S . 
Therefore, from spherical trigonometry 


Sin c - 


sin p _ sin p 

sin (90*- S) cos S 


(4) 


If (x^, y f ) are coordinates in the lunar orbit plane with origin at 
the center of the shadow ellipse, the equation of the shadow ellipse is 


£ + = 


( 5 ) 


-2 



Evaluating Eq (5) at the tangent points. 


X = a. 


— d + R«r sin c 


y'= ± R/r cos c 


gives 


Note that 


From Eq (2) , 


Rj ( I - sin 


+ JSjaincV 

l c)L l a s a s J 


( 6 ) 


-d = 


d -+■ -P 


- d = 


_ -f ~ d _ 


£ [sin ( S - p) 


- , ' . I 

sin (Up)J 


and 


, _ o COsS sin p 

a s -d = — T7 "T" 

sin 5 - Sin P 


- d tan p 

d s £2 n S 


(7) 


Also , 


El 5 Inc - (sin’S- S in* P ) 

Sin S cos £ 

— sine + r — = ^ [l H L-z-(sin*S- sin 1 p')l 

ot s -to^nSL cos X J 


— s*n p cos p 
Sin S cos S' 

( I - sin*c) =• Rf( l_ SS^f) 
Substitution of Eqs (8) and (9) in (6) gives 

2 


( 8 ) 

(9) 


J_ £ 1*2 p \ 

b/ R f n s i n 2 s J 


(10) 


The equation of the ellipse in moon-centered coordinates (x,y) is 


( x + Os — d ) a ,y_ z 

b s a 


a< 


I 


or, in polar coordinates. 


/ r cos cr Q s -d V 

\ a, a s ) 


r* sinV 

b s 2 




The equation for the lunar satellite orbit is 

a ( | - e a ) 


r = 


(ID 


( 12 ) 


I + e COS (cr - CT r ) 
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Then, if values for & and 0}, are specified at any time, a simultaneous 
solution of Eqs (11) and (12) gives the two values of cr at which the satellite 
orbit and the shadow ellipse intersect. Since the equations are transcendental, 
the solution must involve an iterative procedure. Eq (12) may be used to 
eliminate r in Eq (11). 

However, Eq (11) can be simplified without resorting to 
expansions . 


The equation of the form 


(x + 9) : 

a* 


X = | 

b/ 


can be written in polar coordinates as 

(r 2 cosV + 2x3. cos cr ■+ <g a ) r*( 1 - cos*(r) _ 


a. 


b s * 


or 


r z 


cos a cr|| - + r^Zc) cos crj 4 - cj*- a/ = 0 


The discriminate reduces to 


B Z -4AC = 4q/[» + ( 


(13) 


(14) 


From 


d = 


R< 


sin ($ -h p) 

q = a s — d 

and Eqs (3) and (10) , 


9 : 

b 5 a 


b. 2 


co± z £ 


cos*p 


(15) 

(16) 


(17) 


Sin 1 S — sin 2 p 

• _ _^s _ cos 2 S 

b , 2 sin 2 S — sir>*-p 

Substitution of Eqs (15), (16) and (17) into Eq (13) gives the final equation 
of the shadow ellipse in polar, moon-centered coordinates. 

r = R £2$ _ £°s_& sin p ± cos pV I - cos 2 5 cos 1 cr 

COS 2 f=> — cos a S cos a <r 

The positive sign should be chosen in the numerator of Eq (18) . That 
the negative sign is redundant can be seen most easily if the small angle p 
is set equal to zero. Then 


-*( 


V I - cos 1 S 


COS v v 


(p* 


4 - 



Since the radius vector is by definition a positive quantity,, tne 
negative sign is redundant. 


Substitution for r from Eq (12) gives the final equation, to be 
solved for o', the angular coordinate of the intersection of the orbit and the 
shadow ellipse. 

F(<r)= - cos 2 S cos*<r) +[| + 

. \ 

+■ c cos (o'— CoS it cos S sin p -+■ cos pv I - cos*$ cos 2 cr j 

= 0 

This transcesdental equation must be solved by Iteration. Newton's method 
gives 

F( cr n ) 


(19) 


— 


FV„) 


( 20 ) 


where 


F / ( cr )= — Z — — - e - costs’ sin cr coscr 
R a 

— sin ( cr - c r p ) (- cos cr cos S sin p ) 

4- cos — cos 1 $ cos 1 <r ) 4- £ I 

-h e COS ( CT- <7j. sin er oe»s S £ Sin p 
+ cos P cos £ cos cr 1 

y/ I - cos 1 J cosher J 


Equations (19) and (20) then comprise the means of solution for the 
angular positions in orbit at which the orbit intersects the umbra cone. 
The time in shadow can then be determined easily from Kepler's equation, 


X — F. — e. sin E — \ £ (t - t») 

Vo 3 


where 

sin E 
cos E 


\/j ~ e 2 s in O 
| 4 e cos O 

e -+■ cos g 
i e cos G 


6 = cr - cTp 


( 21 ) 


Certain special cases have not been considered in the above solution. 

In Appendix B, the special case in which $<jp and the shadow conic intersection 
is a hyperbola is considered. The solution is shown to be the same a* in the 
elliptical intersection case (Eq 18). Also, in Appendix C the solution for the 
penumbra intersection is shown to be similar to Eq. 18. 
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C. Application to Computation of Lunar Satellite Shadow Time 

The computation routine of Section B was programmed in Fortran IV and 
incorporated in an orbit integration program in order to study shadow 
requirements of the AIMP lunar satellite mission. The integration program 
provides inputs ( S , cr p , a, e as functions of time) to the shadow computation 
routine . 

Orbits to be studied with the shadow routine were generated by the 
following procedure. 

(1) The nominal velocity magnitude, azimuth angle and flight 
path angle at third stage burnout (i.e. injection into the 
earth-to-moon transfer orbit) of the Delta DSV-3E vehicle were 
taken as means of three normal independent distributions of specified 
standard deviations. 

(2) Non-nominal sets of burnout conditions were generated from these 
three distributions for a large number of cases by means of 

a random number generating routine. 

(3) Each set of non-nominal burnout conditions was used as input to 
the Interplanetary Trajectory Encke Method (ITEM) program, which 
for each of the error conditions, as well as initial lunar orbit 
elements resulting from various fourth stage firing times. 

(4) For each transfer orbit one set of initial lunar orbit elements 
which appeared to represent a long lifetime lunar orbit was 
chosen as input to the lunar orbit integration program. If the 
first selection failed to achieve the desired six month lifetime, 
one or two other sets of elements were tested. The lunar orbit 
integration program used is a variation-of -parameters program 
written by Dr. William Kaula with modifications by the Douglas 
Aircraft Company and the Delta Project Office and incorporating the 
shadow routine described in Section B. 

(5) Shadow and lifetime data were tabulated for the lunar orbit having 
the longest lifetime of the cases investigated for each transfei 
orbit. 

Results for a sample of one hundred transfer trajectories are shown in 
the following table. Here, r p , r , a, e, i and co are the initial lunar orbit 
elements (pericynthion and apocynfhion radii in km, semimajor axis in moon 
radii, eccentricity , inclination and argument of pericynthion in degrees 
relative to the lunar equator and vernal equinox) . Also, L is the lunar orbit 
lifetime_in days, t smax is the maximum shadow time experienced in any orbit 
period, is the average shadow time per orbit, and 7. t is the percentage of 
orbit time spent in shadow. The time at which the shadow time first exceeds 
1, 1.5, 2, 2.5 and 3 hours is also given. Here "shadow' 1 means umbra; no 
penumbra investigations were made for this mission. 
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Orbits having apocynthion radii greater than 25,000 nautical miles 
(46,300 km) were not investigated as being of questionable stability. 

The integration program used is not suitable for stability determination 
since accuracy degenerates for high altitude orbits. Also seven transfer 
orbits of the sample did not approach sufficiently dose to the moon to 
achieve a lunar orbit. These cases are indicated in the table. 

Figure 3 summarizes the shadow and lifetime data. For a given time 
after lunar orbit injection (abscissa scale) the ordinate scale indicates the 
number of cases having a longer lifetime and also having a maximum shadow 
time per orbit less than 1, 1.5, 2, 2.5, or 3 hours. The upper curve shows 
lifetime data without maximum shadow time requirements. Figure 4 shows 
histograms of the maximum and average shadow times per orbit. 
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appendix a 


Constants for Cone x.ngle Determination 

Complexity of the shadow time computation can be greatly reduced 
by using mean values of geoid radius, A e , and earth orbit radius, D, in 
the formulas for umbra and penumbra half cone angles (Eqs C.l, C.2). These 
quantities vary between the following limits. 

6256.77 R e 4 6378.165 km 

147.02 x 13 6 4 D 4 152.03 x 10 6 km 

j-'he mean values to he used in general programs may be computed in several 
different ways. 

1. Mean Geoid Radii 


a. Radius of the sphere having volume equivalent to that 
of the oblate snherci^ 


The volume of an oblate spheroid 
(a volume generated by rotating an 
ellipse about its minor axis) is given 
by 


v = 


= 2Tfjy 2 dx 


' f i\ 
V*. j.) 



where : the minor axis coincides with the x-axis. The equation 
of the ellipse in the coordinates shown in the sketch is 


z z 

b 1 + a 2 “ 


(A. 2) 


b. 


Substituting for y in (A.l) gives 

V= lxja'fl- -4)dx = jTa’b 

If R is the radius of the sphere of equivalent volume, 

■j'lrRy 3 = ^Ja x b 

R v - 1 / a* b (A. 3) 

Radius of the ellipse averaged over latitude angle 


In the elliptical cross-section of an oblate 
spheroid, the radius may be expressed as a function 
of latitude as follows: 



- 1 - 


r x cos*(i) r 3 sin 3 4> 

a 1 b 3 “ 


r 2 = 



, si n* d> 
a 1 b s 

The radius averaged over latitude is 

i r* 

R, - ±J r <t)J t • 

(5 

Substitution of (A. 4) in (A. 5) and simplification gives 


(A. j) 


2b 
T j 


\Tl-i 


d 


Then 


where ^ 
-? 


2 cos*<p 


a" 


(A. 6) 


= e 


R * = V'/vrfe; - ^u.f) 


(A. 7) 


where F is the elliptic integral of the first kind, 
For I 


R*~ b[l+2^+9(^) + 50(^-) (A.S) 


c. Radius of the circle having a circumference equivalent to 
that of the spheroid cross-section. 

Determination of the circumference 
of an ellipse is facilitated by 
definition of an auxiliary variable, the 
eccentric anomaly E, and choice of 
coordinates centered at the focus of the 
elliptical cross-section. Then, if s is 
arc length along the circumference, 

ds a = dx 2 -l- dy 2 

x — o. ( cos E — e) dx = — ex si n E d El 

y = b sin E dy = b cos EolE 

ds = V I — e 1 cos 1 El dE 



S = 2a 


j\l I - e 2 cos 2 E dE 

° 'TV 


5= 4 aE (f - 1 ) . J « e 

If the circumference of the equivalent circle is 2 tt R c , 

R c = ^E(f-rj (a. 
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For 


X 


Rc- °['-2i'-3(47-lo(£)-..-] 


(A.IO) 


d. Radius of the circle having an area equivalent to that 
of the spheroid cross-section. 

For the ellipse given by 

L 1 + l l = | 
a 2 + b* 

the area is 


A = I - £ dx 


= -^£x^a a — x 2 -+- cx x sin' 1 U-j 


A = TTab (A. ll) 

If the radius of the circle of equivalent area is R a> 

R« = { ab (a. 12 ) 

Mean Earth Orbit Radius 
a. Astronomical Unit 

The astronomical unit is the mean distance from the sun 
co a fictitious planet, the mass and period of which are those 
used by Gauss in his determination of the solar gravitation 
constant . 


R = 149.53 x 1G 6 km 


Radius averaged over orbit central angle 

r e - || ne)de - 


_ a(l- e 2 ) 


l-e 2 )[ 2 

r L 7T 


2 J. -• \/ 1— e tan 

71 tan 

/l-e* -+• e 


•]; 


r e = a\T7^ 


(a. 13) 


Radius of the circle of area equivalent to that enclosed by the 
earth orbit 


From (A. 11) 

Tab = 7T r a 2 


“ 3 - 



(A. 14) 


r a = fob = c* ( l — e 2 f 

d. Time average of orbit radius 

Kepler's equation is 

1 = E - « sin E 

where 

1 = f£( t - t 0 ) 

Then dE = 1 

di. I - e cos El 

r = a. ( I — e cos ET ) 

dE _ a 

d ’ r 

dE _ nr dE _ i rz 

dt Va 5 d^ " r 1 « 

The orbit radius averaged over time is 

r* = :jr/r( t)dt 

/ o 

where T is the orbit period. 

= - e ; osE) * dE: 
r* = ex ( i + ~ ) (a. 15) 

e. Radius of the circular orbit of equivalent period 

^ . 88 a (A. 16) 

This is also the average obtained by considering the sum of 
the; radii from each focus to every point on the ellipse. 

r, -+- r a = 2 <=» 

£. Radius of the circle of equivalent circumference 
From (A. 10) 
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(A. 17) 



3. Numerical Evaluations 
Mean Geo id Radius 


Type of Average 

Formula 

Value , km 

Sphere of equivalent 
volume 

R v = \/a* b 

6371.02 

Radius averaged over 
latitude 

R.= lpF(i.|) 

6367.45 

Circle of equivalent 
circumference 

R.=*8EU.f) 

6367.47 

Circle of equivalent 
area 

R a — ]j a b 

6367.46 


Assumed constants: 


a = 6378.16 km 
b = 6356.77 km 
e = 0.0818292 


Mean Orbit Radius 


Type of Average 

Formula 

Value, r/a 

Astronomical unit 

Defined, 1.4953xl0®km 


Radius averaged over 
central angle 

r e = ajl- e a 

0.9998601 

Circle of equivalent 
area 

r* = v/Sb = a(l-e a )^ 

0.9999300 

Radius averaged over 
t ime 

r. * a( 1-^) 

1.0001399 

Equivalent period 

r t - a 

1.0 

Equivalent circumferenci 

; r. =^E(i,f) 

0.99993006 


Assumed earth orbit eccentricity e = 0.0167268 
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4. Values for Cone Angle Calculations 

The tables of the previous section indicate that three of the 
four geoid radius averages differ by the order of magnitude of the 
errors in the equatorial and polar radii 

a = 6378.16 + 0.02 km 

b = 6356.77 + 0.05 km 

Therefore based on these values of the physical constants, a suitable 
value for the average geoid radius is 6367.46 km. 

Also, the tables show that all computed values of mean earth 
orbit radii lie within the range of uncertainty in the astronomical 
unit. 

6 6 

a.u. = 149.53 x 10 +0.03 x 10 km, 

which is, therefore, a suitable mean orbit radius. 

For problems requiring greater accuracy than can be achieved with 
mean cone angle values, ephemerides may be used to compute exact 
values for the time in question. Figure A-l shows the variations 
in umbra and penumbra half-cone angles during 1966 for the earth and 
the moon. 
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APPENDIX B 


Special Gas? S<p 


For the case S<p the conic section formed by the intersection of 
the orbit plane and the shadow cone is a hyperbola. However, the equation of 
this shadow hyperbola is the same as that ef the ellipse in the more 
general case £>p(Eq 18). Verification of Eq 18 in the hyperbolic case 
may be carried out in a scheme similar to that of Section B. The geometry 
is shown, in Figure B-l. 


As before, from Figure B-l (A), 

Rt 


d = -: 


sin (p+ S) 


and , from (B) 


sin p 

Sin c = £- 

cos 8 

The equation of the hyperbola in moon-centered coordinates is 


(x - g ') 2 


_y_ 

b a 


(B. 1) 


(B.2) 


(B.3) 


where b 2 = a 2 ( & z — • ) 
c|' — a. •+ d 

The shadow cone is tangent to the me on along a minor circle of latitude 
p above the great circle normal to the M-S line, as in the ellipse case. 
Then any curve lying in. the shadow cone (and, in. particular, the shadow 
hyperbola) is also tangent to the moon on. this circle at 


X, = R „ sine = R ( Sin P 

cos 8 

y i — co s c 

in moon -centered coordinates. Differentiating (B.3) gives 


(B.4) 


At ( V 1» yq) 


2(x-<y) 

a* 

_ I/ dy _ Q 
1? dx U 


b 2 

x-g' 

dx 

a 4 

y 

dy 


CoS C 

dx 


sin P <=> co 3 5 

« . 


(B.5) 


Sin p 


cos c 


sin c - 


i 

Sin C 


— "to 


r\ c. 
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— R„(l •+■ sln a c) -+ 2dsinc 

which in turn, on substitution of (B.l) and (B.2) gives 

a * P? cos p sir»S 
9 sin*p - sin a $ 

Except for sign, this expression is the same as Eq (3) (part B) 

Also, from (B.7) 

I _ q -t- <d — R q sine 
b a a l R< sine 


From (B.9), (B.l) and (B.2) 

a + d — R« sir* C = R<r|— ; 

l Si 

so that (B.10) can be reduced to 


I" sinpcosS sir* p' 
sin 2 p-s'm z S cosSj 


J_ = I /sm*jo _ 

b a R«*lsln*J / 


The equation of the hyperbola may be derived as follows, 
x = r cos c and y - r 8 i n <r , (B.3) gives 


VcosV _ 
\ a* 


— cos ' 1 cr 


Z (c* -4- ol) ^ 

cos cr 


(< 


d 


In the usual quadratic notation, designate 


CO* 


I - C<?S a (T 

b a 


-2(q-fol) 




Since, from (B.9) and (B.ll), 


_L + J _L (si'n a p- sin 2 £) , r 

i_z — c a : — fx r cos b 


«* b * R/ si^^Scos 


so that 


Also, 


Then 


^ __ ^ ( sin p sir* 2 ^ ) / cos* S cosher ^ \ 

Rf s ir> z S \ cos *p I 

B* - 4AC = £[i + sinV(^- d b ^ - - ij] 

r = -B ±fe 2 - 4 AC 
2 A 


(ck + < d) cos cr ± a V 1 "+■ sin 2, cr 

/'2ad-*-d i ,\ 

V b l 7 


.fi*-»-s») co5 * cr - b 1 ] 



From (B.l), (B.9) and (B.10) 
2od + dl 


£-2 | = co-tH 


and 


a -+■ d 


Sin p cos S 


sin ( p ■+■ S) sin (p- S) 

Substitution of (B.14) and (B.15) in (B.13) yields 


o “ sir * P 006 S cos a + COS pV I — cos*£ cosher 
r = rs ^ 

cos*p — cos 4 S c.o* a <r 

This equation is exactly the sans as Eq (18) , Section B. 


(B.13) 

(B.14) 

(B.15) 

(B. 16) 
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APPENDIX C 


Penumbra Determination 


The analysis of Section B is applicable to determination of 
time spent in penumbra as well as time spent in umbra. Figure C-l 
shows the geometry of the penumbra problem. Here (A) shows the penumbra 
and umbra shadow-cones , with their half-cone angles p*, and p, respectively. 


sin p • 



(C.l) 


sin p' „ 


(C.2) 


Values for the earth and the moon are given in the following table 



P> dag 

P 1 , deg 

earth 

moon 

0.2644366 

0.2662145 

0.2693218 

0.2675467 


The geometry of the penumbra problem is obviously the same as that of the 
umbra problem. The magnitude of the half-cone angle changes from p to p , 
and the vertex of the umbra cone is on the opposite side of the central 
body from the vertex of the penumbra cone. In fact, the analysis of 
Section B holds in the penumbra case if the transformation 


P — 



is made. 

This solution can be verified from Figure C-l. 
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The substitution p'-*-—p can be seen to transform (C.l), (C.2), (C.3), (C.4), 
(C,5), (C.6) and (C.7) into the corresponding equations of the umbra problem 
(Section B, Eqs (1),(6)). 


The equation of the shadow ellipse in moon coordinates is 


( X - £ )* 


b/ 


However , 


a-i-f = -a + 2a-d=a 


so that the ellipse equation becomes 


(x + q,-d) * + = 


a. 


b‘ 


5 *->S 

which is the same as the equation preceding Eq (11), Section B. 
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Normalized Circular Orbit Radius (er or mr) 
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